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Abstract 

We study numerically the phonon dispersion relations and corresponding displacement fields for a circular 
cross-section nanowire superlattices consisting of anisotropic GaN and AIN. We determine a set of parameters which 
gives complete phononic bandgaps. The results suggest the potential for manipulating phonons in the micro/nano 
electromechanical systems. 
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Background 

Phononic crystals (PCs) are composite materials made of 
arrays of constituents embedded in host materials [1-3]. 
The interesting characteristics of the PCs are related to the 
existence of phononic bandgaps (i.e., frequency gaps) due 
to the Bragg reflections of the phonons with long wave- 
lengths. We can regard the PC as an opaque barrier for the 
phonons within the phononic bandgaps [4] . This suggests 
the potential for designing various phonon optic devices, 
such as phonon niters, mirrors, resonators, etc. 

Recent advances in fabrication methods enable realiza- 
tion of one-dimensional hetero-structures, i.e., nanowire 
superlattices (NWSLs) [5-10]. Their electronic and opti- 
cal properties were studied, and a variety of possible 
applications utilizing the characteristics were also pro- 
posed [11-14]. In addition, the NWSLs are expected to 
yield interesting physical effects on phonons, which influ- 
ence the electronic states and the transport properties via 
the electron-phonon interaction. These NWSLs can be 
regarded as wire-type phononic crystals ( WPCs), in which 
the phononic bandgaps are induced by the periodicity 
along the wire axis. 

In a previous paper [15], we developed a numerical 
method to derive phonon modes in a free-standing NWSL 
of anisotropic material with an arbitrary shape of cross- 
section. As examples, the phonon modes were calculated 
for the rectangular and square cross-section GaAs/AlAs 
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and InP/InAs NWSLs composed of anisotropic materials 
[15,16]. 

Though above result revealed the important aspects of 
phonon modes in the NWSLs, it seems to be difficult to 
design WPCs with complete phononic bandgaps because 
in the dispersion relations of these NWSLs, many sub- 
bands are folded into the mini-Brillouin zone and the 
frequencies of gaps are different with phonon modes. 
In addition, the gap widths are narrow in these NWSLs 
because the difference of the acoustic impedance between 
the GaAs and AlAs layers (or between the InP and InAs 
layers) is small. The NWSLs with large acoustic mismatch 
would be suitable for designing the phonon optic devices. 

In the present work, we numerically calculate the dis- 
persion relations and corresponding displacement fields 
for a circular cross-section NWSL consisting of GaN and 
AIN, and we determine a set of parameters which gives 
complete frequency gaps. 

Methods 

The equation giving the eigenfrequencies of phonon 
modes in a freestanding NWSL composed of anisotropic 
crystals was formulated in [15]. In this method, the dis- 
placement components Ui (i = x,y,z) are expanded in 
terms of a set of basis functions (p a (r) 



Ui(r) = ^Wa(r). 



(1) 
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The expansion coefficients A a { and the eigenfrequencies 
co are determined by solving the generalized eigenvalue 
equation: 



(2) 



For the NWSLs composed of cubic materials, the matrix 
elements Hpi&i and S^ a t are written as 



Here, C^ v and p are the stiffness tensor and mass den- 
sity, respectively, which are dependent on r in the NWSLs. 

As basis functions, we adopt the product of powers of 
the Cartesian coordinates in the xy plane. The z depen- 
dence is expressed in the form of the Bloch wave: 



4>a (k,r) 



1 /xy* 



(k+G)z 



(7) 



Hpuat = | W \C n I to) + W l c 44l/ a ) 



&u + m\c 12 \tcc) 



+ {ip\C 44 \ia)}(l-8 it ), 



Spi,at =8it(P\p\a), 
where 



mC»Aka) = j^C^-^dr, (5) 



(P\P\«) 



4>Ur)p(r)(j) a (r)dr. 



(3) 



(4) 



(6) 



Here, R denotes the radius of the wire; G is the recip- 
rocal lattice vector; and V = ttR 2 D is the volume of the 
unit cell, where D is the length of the unit cell in the 
z direction. The basis functions are specified with a = 
(m, n, G). The expressions for the matrix elements Hp^i 
and Spi }0l i can be analytically obtained for the circular 
NWSLs.' 

Based on group theory [17,18], the phonon modes are 
classified symmetrically. In the present study, the con- 
stituent layers are assumed to be cubic materials, i.e., 
zinc-blende structure. The group of k is for 0 < k < 
7i ID. The irreducible representations are Ai, A2, Bi, B2, 
and E [15]. 

By considering the above symmetry, the symmetry- 
adapted basis function in the present system can be con- 
structed, and the phonon dispersion relations of each 
mode are independently calculated. 




k 



Figure 1 Phonon dispersion relations of the circular cross-section GaN/AIN NWSL. (a) R = 5.0 nm and D = 20.0 nm; (b) R = 5.0 nm and D = 45.0 
nm (thicknesses of GaN and AIN layers are 25.0 and 20.0 nm, respectively). The B2 modes exist in higher frequency range, and all the E modes are 
doubly degenerated. 
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Results and discussion 

The boundary condition at the free surface of the wire 
requires that the wave numbers in the lateral direction 
are discretized. On the other hand, the wave num- 
ber k in the longitudinal direction has a continuous 
value. Therefore, even the homogeneous plain nanowire 
has subband structure. In the NWSL, the subbands are 
folded into the mini-Brillouin zone. In other words, 
the size of the mini-Brillouin zone and corresponding 



phononic bandgaps are determined by the periodicity D 
of the NWSL. 

For the nanowires with smaller R, the maximum wave- 
lengths in the lateral directions become shorter. As a 
result, the subbands except for the lowest dispersion curve 
of each mode go up to the higher frequency region, though 
a lot of subbands are folded into the mini-Brillouin zone. 
Changing the ratio of R and D, we can control the phonon 
modes in the lower frequency range. 



(a) Ai mode (b) A2 mode 





Figure 2 Displacement field patterns corresponding to the lowest (a) Ai , (b) A2, (c) Bi , (d) B2/ and (e) E modes at k = it ID. 
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Figure la illustrates the phonon dispersion relations 
calculated for R = 5.0 nm and D = 20.0 nm (thicknesses 
of both GaN and A1N layers are 10.0 nm). Other param- 
eters we used are as follows: p = 6.15 g/cm 3 , Cn = 
2.96, C12 = 1.54, and C 44 = 2.06 (all in units of 10 12 
dyn/cm 2 ) for GaN; p = 3.26 g/cm 3 , Cn = 3.04, C 12 = 
1.52, and C44 = 1-99 (all in units of 10 12 dyn/cm 2 ) for 
A1N [19]. 

In the present frequency range, we can see four dif- 
ferent modes, i.e., the Ai, A2, Bi, and E modes. The 
dispersion curves corresponding to the B2 modes exist 
in higher frequency range. The lowest dispersion curve 
of the A2 mode is mostly overlapped with that of the 
Ai mode. 

The complete frequency gaps are realized for the param- 
eters we selected. For comparison, we show in Figure lb 
the phonon dispersion relations calculated for R = 5.0 nm 
and D = 45.0 nm (thicknesses of GaN and A1N layers 
are 25.0 and 20.0 nm, respectively). In this example, the 
broader bandgaps of the Ai, A2, and E modes are nearly 
coincident with each other. Here, we note that the com- 
plete bandgaps disappear if the isotropic approximation is 
used for each constituent layer. 

The lowest dispersion curves of the Ai and A2 modes 
are linear in /<, i.e., co vanishes at k = 0. On the other hand, 
the lowest frequencies of the Bi and B2 modes have finite 
values. For thin (thick) NWSLs, the lowest frequencies of 
Bi and B2 modes become higher (lower). 

Figure 2a,b shows the displacement patterns corre- 
sponding to the lowest Ai and A2 modes at k = tt/D, 
respectively. The lowest Ai mode has the feature of a 
dilatational mode, while the lowest A2 mode shows the 
feature of a torsional mode. 

Figure 2c,d,e shows the displacement patterns corre- 
sponding to the Bi, B2, and E modes at k = tt/D, 
respectively. The lowest Bi mode shows the feature of a 
stretching mode, i.e., the alternating dilatation and con- 
traction in the x and y directions, while the lowest B2 
mode shows the features of a shear mode, i.e., alternating 
stretching in the two diagonal directions. 

All dispersion curves of the E mode are doubly degener- 
ate because the irreducible representation of the E mode 
is two-dimensional. For the E modes, the lowest disper- 
sion curve near k - 0 is proportional to k 2 . This parabolic 
behavior is due to the fact that these modes correspond to 
the bending of the NWSL. Figure 2e clearly shows that the 
E modes have a feature of flexural mode. 

Conclusion 

We theoretically studied the acoustic phonon modes in 
circular cross-section NWSLs consisting of cubic GaN 
and A1N. We calculated their dispersion relations and 
phonon displacement fields. These modes are classified 
into five types, i.e., Ai, A2, Bi, B2, and E modes, which 



have features of dilatational, torsional, stretching, shear, 
and flexural modes, respectively. We determined a set 
of parameters which gives complete phononic bandgaps. 
The results suggest the realization of the optimized 
phonon devices, such as phonon filters or mirrors in the 
micro/nano electromechanical systems. 

In the present work, we only showed the results for 
the NWSLs consisting of cubic materials (i.e., zinc-blende 
structure). The results for the NWSLs with wurtzite struc- 
ture will be given and the difference will be discussed 
elsewhere. 
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